Strong enhancement of the d-wave superconducting state in the three-band Hubbard 

model coupled to an apical oxygen phonon 
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We study the hole binding energy and pairing correlations in the three-band Hubbard model 
coupled to an apical oxygen phonon, by exact diagonalization and constrained-path Monte Carlo 
simulations. In the physically relevant charge-transfer regime, we find that the hole binding energy 
is strongly enhanced by the electron-phonon interaction, which is due to a novel potential-energy- 
driven pairing mechanism involving reduction of both electronic potential energy and phonon related 
energy. The enhancement of hole binding energy, in combination with a phonon-induced increase 
of quasiparticle weight, leads to a dramatic enhancement of the long-range part of d-wave pairing 
correlations. Our results indicate that the apical oxygen phonon plays a significant role in the 
superconductivity of high-Tc cuprates. 



I 

O 

o 



> 
O 

(N 

^' 

O 
O 



X 



I. INTRODUCTION 

Despite years of intensive research, the pairing mecha- 
nism responsible for d-wave superconductivity (dSC) in 
the high-Tc cuprates remains a puzzle. It is generally 
believed that the conventional phonon mechanism is in- 
consistent with d-wave pairing symmetry and not strong 
enough to explain transition temperatures higher than 
100 K. Most investigations in this direction have been 
focused on the pure electronic mechanism, but no con- 
sensus has been reached so far. Recently, accurate exper- 
iments displayed pronounced phonon and electron-lattice 
effects in these materials, which are manifested by a large 
softening and broadening of certain phonon modes in the 
whole doping region. In particular, the in-plane copper- 
oxygen bond-stretching phonon, apical oxygen phonon 
(AOP), and oxygen Big buckling phonon are shown to 
be strongly coupled to charge carriersii^. Moreover, 
photoemission-spectroscopy-resolved kink structures — 
are probably caused by coupling of quasiparticles to 
phonon modes. These findings suggest that phonons are 
important for the physical properties of high-Tc cuprates. 

Various theoretical attempts have been made to un- 
derstand the role of phonons in high-Tc superconductiv- 
ity (HTSC)^"— , but the answer remains unclear. In a 
functional renormalization group study^, Honerkamp et 
al. found that the Big buckling phonon enhances the 
d-wave pairing instability in the Hubbard model. More 
recently, an exact diagonalization (ED) study of the t — J 
model coupled to phonons shows that coupling to the 
buckling mode stabilizes d-wave pairing while coupling 
to the breathing mode favors a p-wave pairing^. On 
the other hand, based on dynamical cluster Monte Carlo 
calculations of the Hubbard model coupled to Holstein, 
buckling and breathing phonons, Macridin et al. found 
that while these phonons can indeed enhance pairing, a 
strong phonon-induced reduction of quasiparticle weight 
leads to a suppression of dSCS. 

In this paper, we study the effect of AOP on dSC in 



the more realistic three-band Hubbard model. Our work 
is motivated by recent experiments showing that the dis- 
tance between apical oxygen and the Cu02 planed, the 
disorder around apical oxygeniSiii, and the apical hole 
stated have significant effects on Tc. Basically, there 
are two routes for apical oxygen to affect HTSC: one is 
to tune the electronic structure of Cu02 plane, leading 
to a change of T^; the other one is to directly couple 
apical oxygen vibrations to charge carriers on the con- 
ducting Cu02 planei^ii^. Here, we focus on the second 
route and study a strongly anharmonic vibration of apical 
oxygen in a double-well potential, which is evidenced in 
the X-ray absorption spectroscopy of several typical high- 
Tc compoundsi^^— . For a strongly anharmonic motion 
in the double-well potential, the first excitation energy 
Ai? = El — Eq is much smaller than the ones excited to 
higher energy levels En>2- In this case one can take into 
account only the lowest quantum states $o and <f>i with 
energies Eq and Ei and model the low-energy motion of 
apical oxygen by a local two-level system represented by 
a pseudospi n^^i^" degree of freedom. 

Our main results, obtained by ED and constrained- 
path Monte Carlo (CPMC) methods, arc presented in 
Figs. \V[h), (Ha) andU Fig. [i;a) clearly shows that the 
coupling to the AOP induces a strong enhancement of 
hole binding energy, and this enhancement effect grows as 
the Coulomb repulsion Ud on the copper site is increased. 
An analysis of the contribution of different energies to 
the hole binding energy reveals a novel potential-energy- 
driven pairing mechanism that involves reduction of both 
electronic potential energy and phonon related energy. 
As a combination of increasing pairing interaction and 
quasiparticle weight (sec Fig. (Ub)), the d-wave pairing 
correlations arc found to be strongly enhanced by the 
electron-phonon (cl-ph) coupling (see Fig.|3|). 

Our paper is organized as follows: In Section |TT1 we 
define the Hamiltonian and the physical quantities cal- 
culated and discuss the choice of model parameters. In 
Section [nil we present our numerical results and discuss 
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the physical mechanism responsible for the AOP-induced 
enhancement of dSC. Finally, in Section HVl we discuss 
in detail our main conclusions. 



II. MODEL AND NUMERICAL APPROACH 

To model the electronic structure of Cu02 plane and 
the coupling of holes to the anharmonic AOP, we adopt 
the following Hamiltonian proposed in Ref.— , 



H — Hk + Hpot + Hph , 



(1) 



where Hk, Hpot, and Hph stand for the kinetic motion of 
holes, the potential energy for holes, and phonon related 
energy, respectively. They are expressed in the form: 



+ tipip]aPka+h.C.), 
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Hpot = eE^^'o + ^-^E 
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(2) 



(3) 



and 



against ED on the 2x2 cluster, and the difference for 
the electronic kinetic energy, as well as for the charge and 
magnetic moment at the copper sites, is less than 3% up 
to Ud = 6. 

The hole binding energy is defined as: 



A = E2 + Eo-2Ei, 



(5) 



with En the ground-state energy for n doped holes. The 
(ia;2_j,2 pairing correlation is defined by. 



Pd{R) = (At(i?)A,(0)), 



(6) 



where 



^ ^Rt R+Si ^Ri^R+St" 



with 6 = ±x, ±y. fd{S) = 1 for (5 = ±i and fd{S) = — 1 
for 5 = zty. We calculate also the vertex contribution to 
the correlations defined as follows: 



Vd{R) = PdiR) - PdiR), 



(7) 



where Pd{R) is the bubble contribution obtained with 
the dressed (interacting) propagator—. 



Here, the operator creates a hole at a Cu 3dj.2_y2 
orbital and creates a hole in an 02px or 2py orbital. 
sf and sf are the pseudospin operators for S* = 1/2. Ud 
denotes the Coulomb energy at the Cu sites, t^^ = ±tpd 

and t^pp = ±tpp are the Cu-0 and 0-0 hybridizations, 
respectively, with the Cu and O orbital phase factors in- 
cluded in the sign. The charge-transfer energy is e, i.e. 
the oxygen orbital energy, gcu and go in Hph denote the 
strength of el-ph coupling. 51 stands for the tunneling 
frequency of the two-level system. S' is the vector con- 
necting Cu and its nearest-neighbor (NN) O. According 
to quantum cluster calculations^, the parameters lie in 
the range: tpd = 1.3 — 1.5eF, tpp = O.GeV, e = 3.6eV, and 
Ud = 8.5 — 10.5eV. In units of tpd, we choose a parameter 
set tpp = 0.3 and e = 3, while Ud is varied from weak to 
strong coupling, including the physical value Ud = 6 — S- 
51 = 0.5 and gcu ^ go = g are assumed for the results 
presented below, except explicitly noted otherwise. 

Our calculations arc performed on clusters of 2 x 2, 
8 X 4, 6 X 6 and 8x6 unit cells with periodic bound- 
ary conditions using the ED and CPMC method o^^'^'^ . 
In the CPMC method, we follow Refs. [H and [ll to use 
the Worldline representation for the pseudospins and pro- 
jected the ground state j^E'o) of el-ph interacting system 
from a trial wave function |4't) = I^J) ® I^D- Here, 
l^*^) and represent the hole and phonon parts, 

respectively. The CPMC algorithm has been checked 
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RESULTS AND DISCUSSIONS 



First, we show ED results for the 2x2 cluster with 
one hole doped beyond half filling. The electronic ki- 
netic energy Ek = (Hk), the peak value A''"^"-'^ of the 
single-particle spectral function ^1(^,0;) at the Fermi 
energy (w = 0), the charge ricu ~ (nf) and the 
magnetic moment = {{nf^ — nf^)'^) at the copper 
sites are displayed in Figs. [Ua)-md). Here, A{k,uj) = 

-7 E™('fo(l)|c™,fe.(w-i5i+ff+z77)-ic|„^,J$o(l)) with 
$0(1) and El denoting the ground-state wave function 
and its energy for one-hole doping. The index m corre- 
sponds either to the (1^2 _y2 or to the Px,y orbitals, and 
?7 = 0.2. One can clearly see that Ek is lowered with 
increasing the el-ph coupling g at all Coulomb energies. 
Meanwhile, an increase of A^'^"'^ with increasing g indi- 
cates that the quasiparticle weight is increased by the 
el-ph coupling. In contrast, previous studies of harmonic 
phonons in the Holstein-Hubbard model found that the 
kinetic energy of electrons is increased with increasing el- 
ph coupling, accompanying a reduction of quasiparticle 
weight^i2^. To explore the physical reasons for lowering 
Ek, we switch off the coupling of apical phonon either 
to copper or to in- plane oxygen, i.e., we set g^u = or 
(7o = in Eq.(|3]). It is found that Ek is lowered for the for- 
mer case, but increased for the latter case. These results 
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FIG. 1: (color online)(a) Electronic kinetic energy Eh (per 
unit cell); (b) peak value JiJ"^°-^ of the spectral function 
A{k = {■K,Q),uj) at the Fermi energy; (c) charge ricu and (d) 
magnetic moment at the copper sites as a function of Ud 
at different el-ph coupling g. The inset of figure (b) displays 
A(k = (7r,0),aj) as a function of uj at Ud = 6, with g = 0.0 
(solid line) and g — 0.8 (dashed line). The The results are 
obtained from ED for the 2x2 cluster. 



demonstrate that the special coupling of apical oxygen 
phonon to in-plane oxygen is responsible for lowering the 
electronic kinetic energy. 

From Fig.[ljc), we notice that the charge is transferred 
from copper to oxygen sites, which, in combination with a 
phonon-mediated retarded attraction between holes with 
opposite spins, results in a reduction of magnetic moment 
at Cu sites, as shown in Fig.[ljd). Similar effects of AOP 
on Eti, ricu and are also observed for larger clusters 
obtained by CPMC simulations, and representative re- 
sults on the 6x6 cluster are shown in Table ID The last 
column in Tab. |I] shows that the value of NN Cu — Cu 
spin correlation (Sf • 5*1) becomes less negative with in- 
creasing g, implying a suppression of antiferromagnetic 
(AFM) spin correlation. 



The hole binding energy A is shown in Fig. [2ja) as a 



TABLE I: g dependence of Ek (per unit cell), ricu, rn? , and 
NN Cu — Cu spin correlation {Si ■ Sj) on the 6x6 cluster 
at C/d = 6. The number of holes Nh = 42, corresponding to 
a hole doping density x ~ 0.f67. Statistical errors are in the 
last digit and shown in the parentheses. 
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0.0 


-2.612(2) 


0.7563(1) 


0.6937(1) 


-0.0800(3) 
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-2.659(3) 


0.7426(2) 


0.6786(1) 


-0.0715(3) 


0.4 


-2.691(6) 


0.7292(3) 


0.6651(2) 


-0.0659(6) 


0.6 


-2.723(9) 


0.7142(5) 


0.6505(5) 


-0.0605(9) 



FIG. 2: (color online) (a) Hole binding energy A obtained 
from ED for the 2x2 cluster as a function of Ud- (b)-(d) show 
the contributions to A from different energies, as discussed in 
text. The el-ph coupling g is indicated by the shape of the 
symbol. 



function of the Coulomb energy Ud at different g. At all 
Ud, the binding energy is decreased by switching on the 
el-ph coupling, signaling an enhancement of hole pair- 
ing interaction. It is remarkable that this enhancement 
effect becomes stronger with increasing Ud, which is par- 
ticularly evident in the region Ud < S- 

In order to identify the physical origin for this enhance- 
ment, the contributions to A from the hole kinetic en- 
ergy Ek, the hole potential energy Epot = {Hpot), and 
the phonon related energy Eph = {Hph) are depicted in 
Figs. [2] (b)-(d), respectively. Here, A^, Apot and Aph 



have similar definitions to A, with E in Eq.(l5l) replaced 
with Ek, Epot and Eph, respectively. These quantities 
represent the gain in the corresponding energy when the 
second hole is doped in the vicinity of the first one. Al- 
though the kinetic energy of holes is reduced upon in- 
creasing the el-ph coupling (see Fig. IHa)), an increase 
of Afc with increasing g displayed in Fig. [2][b) indicates 
that the kinetic energy gain for two doped holes is re- 
duced by the el-ph coupling. As seen in Fig. [2jc) and 
Fig. EJd), a decrease of Apot and Aph with increasing g 
reveals that it is the reduction of electronic potential en- 
ergy and phonon related energy between two doped holes 
that enhances the hole binding energy A. In addition, the 
decrease of Apot and Aph becomes more pronounced as 
Ud is increased, leading to a stronger enhancement of the 
hole binding energy in the strong correlation regime (see 
Fig. [Ha)). 



Fig. [3] displays the hole binding energy and different 
contributions to A on the 2x2 cluster with antiperiodic 
and mixed (periodic in the x direction and antiperiodic 
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FIG. 3: (color online) (a) Hole binding energy A obtained 
from ED for the 2x2 cluster as a function of Ud- (b)-(d) show 
the contributions to A from different energies. Open circle is 
for g = 0, and filled symbols for g = 0.4 with antiperiodic 
(square) and mixed (diamond) boundary conditions, respec- 
tively. 



in the y direction) boundary conditions. A comparison 
of the results in Figs. [2] and [3] shows that although the 
boundary condition has strong effects on the amplitude 
of hole binding energy, the AOP-induced enhancement of 
hole binding energy is qualitatively similar for different 
boundary conditions. This demonstrates that our find- 
ings reflect the intrinsic effects of AOP in the studied 
model. 



Based on the ED results on the small cluster, we can 
conclude that the coupling of AOP to holes can enhance 
superconductivity on the Cu02 plane. The question aris- 
ing is whether the d-wave pairing symmetry is enhanced? 
This issue can be addressed by examining the behavior 
of the d-wave pairing correlation Pd{R) in Eq.([ni). In 
Figs. [Ha)-(c) wc show Pd{R) as a function of R for the 
6x6 cluster at Ud = 2, 4 and 6, respectively. The 
corresponding vertex contribution is also displayed in 
Figs. [l{d)-(f). In the weak-correlation case {Ud = 2), 
we observe that Pd{R) is modified slightly by the el-ph 
coupling. However, when Ud is increased to 4 and 6, the 
pairing correlation is enhanced at all long-range distances 
for R > 2. At g — 0.6, the average enhancement of the 
long-range part of Pd is estimated to be about 52% and 
65% for Ud = 4 and Ud = Q, respectively. This increasing 
enhancement of Pd with Ud is consistent with our find- 
ings for the binding energy and quasiparticle weight. A 
dramatic increment of the vertex contribution with in- 
creasing g, as shown in Fig. HJ^e) and Fig. IHf), further 
provides strong evidence that the d-wave pairing interac- 
tion is actually increased by the cl-ph coupling. 
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FIG. 4; (color online) Pd{R) as a function of the distance R 
between hole pairs on the 6x6 cluster at Ud = 2 (a), Ud ~ 4: 
(b) and Ud = 6 (c). (d)-(f) show the vertex contribution 
Vd{R) corresponding to the left panels. Results are obtained 
for Nh = 42 holes. The el-ph coupling g is indicated by 
the shape of the symbol. The inset of figure (c) displays the 
average Pave (see text) as a function of Q. 



We also study the effect of the tunneling frequency 
on the pairing correlations. In the inset of Fig. |4] (c), the 
average of long-range d-wave pairing correlation. Pave ~ 
{1/N') J2r>2 Pd{R) where N' is the number of hole pairs 
with i? > 2, is plotted as a function of Vt for Ud = & and 
g = 0.4. Wc notice that its frequency dependence is 
rather weak, suggesting that the isotope effect of apical 
oxygen on the superconductivity is very small. This is 
in good agreement with the site-selected oxygen isotope 
effect in YBa2Cu30e+J^. 



To illustrate the effect of hole doping on the plionon- 
induced enhancement of dSC, in Fig. [5] we show Pave as 
a function of hole doping density x at Ud = 6 for the 
6 X 6, 8 X 4, and 8x6 clusters. The combination of the 
results on the three clusters shows that in a wide hole 
doping region 0.0 < x < 0.5, the d-wavc pairing correla- 
tion is enhanced by the el-ph coupling, and as the hole 
doping density is increased from underdoping to optimal 
doping and then to overdoping, the enhancement of Pave 
is weakened monotonically. The reduced enhancement of 
Pave with larger hole doping, together with the stronger 
enhancement of A and Pd with increasing Ud, demon- 
strates that strong electronic correlations and/or AFM 
fluctuations play a crucial role in the phonon-induced 
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enhancement of dSC. 
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FIG. 5: (color online) Pave vs the hole doping density x for 
the 6 X 6 (a), 8 X 4 (b), and 8 x 6 (c) clusters at Ud = 6. 



IV. 



CONCLUSIONS 



In summary, our numerical simulations show that the 
hole binding energy is strongly enhanced by an AOP- 
induced reduction of electronic potential energy and 
phonon related energy. As a combination of two con- 
curring effects, i.e. the enhancement of hole pairing in- 
teraction and the increase of quasiparticle weight, the 
long-range part of d-wave pairing correlations is dra- 
matically enhanced with increasing the el-ph coupling 
strength. Our results also show that strong electronic 
correlations and/or AFM fluctuations are crucial for this 
phonon-induced enhancement effect. The consistent be- 
havior of our results on different clusters suggests that 
the phonon-induced enhancement of dSC could survive 
in the thermodynamic limit. 
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